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In this paper we determine the exact fermionic spectral function of the Bloch-Nordsieck model at 
finite temperature. Analytic results are presented for some special parameters, for other values we 
have numerical results. The spectral function is finite and normalizable for any nonzero temperature 
values. The real time dependence of the retarded Green's function is power-like for small times and 
exhibits exponential damping for large times. Treating the temperature as an infrared regulator, 
'^^^ ^Iso give a safe interpretation of the zero temperature result. 

^ ' 

o : 

(N ■ 

C3 ! I. INTRODUCTION 



0^ 



The behavior of the ultra-soft regime of massless field theories presents a serious challenge which, on the other hand, 
is crucial for understanding of the most, physically relevant theories. The soft nature of the excitations perturbatively 
— I , leads to infrared (IR) divergences in various physical quantities such as the self-energy near the mass shell. To have 
^ ' reliable results, one has to resum the most sensitive part of the IR physics. The identification of the sources of these 
, divergences, the elaboration of the appropriate mathematical tools and finally the realisation of the resummation itself 
is a formidable task. Moreover, the details can depend on the environment, that is the resummation in deep inelastic 
^ \ scattering and at finite temperature equilibrium may require diff'erent approaches. It is not a surprise, therefore, that 
so many resummation methods exist, working at different circumstances. 

In this sensitive field, where the physical reliability of a resummation may crucially depend on the correct identifica- 
tion of the relevant sources of the IR divergences, it is quite valuable to find a model which is physically motivated and 
exactly solvable. This is the reason why so many two-dimensional conformal and integrable theories have relevance. 
, In four dimensions, however, exactly solvable models are much rarer. 

A physicallv well- motivated model is the Bloch-Nordsieck (BN) model [1] . In its long history it became a text-book 
material [H, yj. Physically it corresponds to the deep IR limit of QED, where the photons have no energy even for a 
fermion spin-fiip. In particular it can be used to prove QED theorems in this energy regime 3- The BN model can 
be solved exactly, the photon contributions can be fuUv summed up. The fermion propagator has been calculated at 
zero temperature both based on functional methods ^^'i with help of Dyson- Schwinger equations [1, 0]; where 

also a detailed renormalization analysis is possible. The spectral function of the model at zero temperature reads in 
Feynman gauge as 
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q{'w, T ^ 0) ^ w ^ ^ , where w — u^^-p^ — m. a = — . (1) 
H ■ ' 47r 



Here m is a 4-vector parameter of the model, loosely identifiable with the four- velocity of the fermion. This function, 
however, has a singular behavior: it is not normalizable, therefore the sum rule J p — 2t: can be satisfied only 
with zero wave function renormalization factor. Moreover, the naive inverse Fourier transform of this function is 
~ t"/'^ describing growth of correlation in time. The correct physical interpretation of these results requires some IR 
regulator, which can be, for example, the temperature. 

At finite temperature the model is much less studied. In the seminal papers of Blaizot and lancu 0, Q the authors 
studied the large time behavior of the fermion propagator with the Hard Thermal Loop (HTL) improved photon 
ropagator. Using this result, Weldon worked out a spectral function which is valid in the vicinity of the mass shell 
With a different approach. Fried et. al. studied the time dependence of the momentum loss of a hard incoming 
fermionic particle |1C|. 

We have several goals in this paper. The main goal is to work out the complete spectral function of the BN model 
for all momenta, and see how the short time dynamics, resembling the T — ^ limit, goes over to the long time 
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damping. Because of the relative simplicity of the model we can even give analytic solutions for certain parameters, 
while for other, analytically not reachable parameter values we used a well controlled numerical procedure. Another 
goal is to extend our Dyson-Schwinger formalism combined with Ward identities , which works excellently at zero 
temperature, to finite temperatures. With the help of it, the complete renormalization process remains fully controlled. 

Our paper will be organized as follows. First, we define the Bloch-Nordsieck model in Section |lll We review the 
Dyson-Schwinger equations and Ward identities in finite temperature real time formalism, and apply them to the 
Bloch-Nordsieck model. In Section IIIII we solve these equations. At zero velocity (Subsection IIII C[) we provide an 
analytic formula for the fermion propagator, supported by a numerical verification. At nonzero velocity (Subsection 
IIIID|) we solve them numerically. In Subsection IIII El we compare our results with previous works in the literature. 
In Section HV] we give the conclusions of the paper. 



II. THE BLOCK NORDSIECK MODEL AT FINITE TEMPERATURE 



The Bloch-Nordsieck model is the low energy limit of QED, where we take into account only a single spin orientation. 
Its Lagrangian is related to the QED Lagrangian by changing the Dirac matrices for a four- vector u^: 

C = -^F^,F^'' + ^\iu^D^' -m)^, iD^ = id^-eA^, F^, ^ d^,A, - d,A^. (2) 

We can choose u to be a four- velocity, or it can be u = (1, v): the two are related by a simple field and mass rescaling, 
since by replacing ^I' ^! / ^Juq and m — muo, we can reach the — 1 scenario. The quantity v = u/uq can be 
interpreted as the velocity of the fermion. 

We are interested in the finite temperature fermion propagator. To determine it, we use the real time formalism 
(for details, see Here the time variable runs over a contour containing forward and backward running sections 

(Ci and C2). The propagators are subject to boundary conditions which can be expressed as the KMS (Kubo-Martin- 
Schwinger) relations. The physical time can be expressed through the contour time t = T{t). This makes possible 
to work with fields living on a definite branch of the contour, 4'a(i,x) — 4'(Ta,x) where T{Ta) — t, and Tq G Ca for 
a = 1, 2; and similarly for the gauge fields. The propagators are matrices in this notation: 

iGabix) = (Tc^a{x)^l{0)) and zG^.,a6(x) = {TcA^a{x)A,b{0)) , (3) 

where Tc denotes ordering with respect to the contour variable (contour time ordering). Gn corresponds to the 
Feynman propagator, and, since the C2 contour times are always larger than the Ci contour times, G21 = G^ and 
6*12 = are the Wightman functions. The KMS relation for a bosonic/fermionic propagator reads Gi2(i,x) = 
±G2i(t — which has the following solution in Fourier space 

iGi2(fc) = ±n±{ko)Q{k), iG2i{k) = (1 ± n±)(/co)f?(/c), (4) 

where 

n±{ko) ^ — - and g{k) ^ iG2i{k) ~ iGi2{k) (5) 

are the distribution functions (Bose-Einstein (-I-) and Fermi-Dirac (-) statistics), and the spectral function, respectively. 
It is sometimes advantageous to change to the R/A formalism with field assignment ^'1^2 = ^'r ± ^'a/2. Then one has 
Gaa = for both the fermion and the photon propagators. The relation between the 12 and the R/A propagators 
reads 

Grr = ^ , Gil = Gra + G12, Q — iGra — iGar- (6) 

The Gra propagator is the retarded, the Gar is the advanced propagator, Grr is usually called the Keldysh propagator. 
At zero temperature the fermionic Feynman-propagator reads: 

Q^{p) = TT^-—-- (7) 

It has a single pole which means that there is no antiparticles in the model. Consequently, all closed fermion loops 
are zero, thus there is no self-energy correction to the photon propagator at zero temperature. Physically this means 
that the energy is not enough to excite the antiparticles. In fact, if we interpret the u parameter as the four-velocity 
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of the fermion, the Bloch-Nordsieck model describes that regime where the soft photon fields do not have energy even 
for changing the velocity of the fermion (no fermion recoil). This leads to the interpretation that the fermion is a 
hard probe of the soft photon fields, and as such it is not part of the thermal medium [8]. So we will set Q12 = 0, 
therefore the closed fermion loops as well as the photon self energy remain zero even at finite temperature. Another, 
mathematical reason, why we must not consider dynamical fermions - which could show up in fermion loops - is that 
the spin-statistics theorem [l3] forbids a one-component dynamical fermion field. 
This means that now the exact photon propagator reads in Feynman gauge 



, ^?(fc) = 2^sgn(fco)<5(fc2), (8) 



all other propagators can be expressed using identities (|4|) and dH). 

A. Dyson- Schwinger equations 

The operator equations of motion give relations of the different Green's functions, formulated as the Dyson- 
Schwinger equations. These equations are locals and so they are valid in generic non-equilibrium situations, and, 
of course, in a thermal medium, too. 

The generating form of the Dyson-Schwinger equations for generic fields reads [l2j | 

/ SS \ " 



In real time formalism the time variable is the contour time (usually it is the variable of the path integral). We define 
the fermionic self energy in the usual way 



g{x, y) = {x, y) + / ct'x'dy G^'^ (x, x'^x' , y')g{y', y), (10) 
Jc 

where the symbol J^^ means time integration over the contour. Then we find in the Bloch-Nordsieck model 

'E{x,y) — ia{xo)e'^Ufj_ / d'^wd'^zQ{x,w)G'^'^ (x, z)Ti,{z;w,y), (11) 
c 



where the tree level vertex is eu^, the proper vertex is denoted by eF^, and a{xo) is 1 if xq € Ci and —1 if a;o £ C2. 
This factor appears because we expressed the functional derivative g^^^y^ through the derivatives of the Lagrangian, 
which, however, changes sign on C2. 

We can also express this equation with the two-component notation as it can be seen on Fig. [TJ In terms of analytic 




FIG. 1. The Dyson-Schwinger equations in real time formalism 

formulas it reads: 

T,ab{x,y) = iaac'^Uf^ d'^wd'^zgac{x,w)G'^'^{x, z)T^.dcb{z;w,y), (12) 

c,d=l 

where aa = In Fourier space it reads: 

^ab{p) = iaas^u^ ^ / — — ^ Gacip ~ k)G^^''^{k)T ^-dcbik] p - k,p). (13) 



B. The vertex function in the Bloch-Nordsieck model 

The second use of the Dyson-Schwinger equation is to have a form for the vertex function. From ^ we find for 
any gauge theories 

where O is any local operator containing ^ and ^. This implies, in particular 

{A^{x)^{y)^{z))^ [ d*x'G,,{x,x'){r{x')^{y)^{z)), (15) 
Jc 

where is the conserved current. The vertex function shows up in the A^t'^E'^ correlator as 

{A^{x)'f{y)^{z)) - / d^x'd^'d^z' iG^,{x, x')ig{y, y'){~ie)V%x' , y' , z')ig{z' , z). (16) 
Jc 

From here we find 

dyd'^z' ig{y, u) er''(x; w, v)ig{v, z) = (j'^(a;)^'(y)§(z)) . (17) 



In the BN model the fermion propagator is a scalar, moreover = eu^^E'^^ is proportional to u^. Therefore the 
vertex function is proportional to u^, too. This is written in the Fourier space as 

r^{k;p, q) = uf^T{k;p, q) {2nf5(k+p- q), (18) 

where we also used the energy-momentum conservation. 

C. Ward identities 

The local equations expressing current conservation can be used in a similar manner. The generating form reads 
d " 

g-^ (i^(a:)$ai (a;i) . . . $a„ (Xn)) ^-^Yl ^^<^Ji^ " ^k) ($ai {Xl) . . . $a._i (Xfe-i ) (y) (Xk+l) . . . $a„ (Xn)) , 

^ k=l 

(19) 

where A$i is the transformation of the ith field generated by the conserved charge Q — J d^xj'^{t,x.). This means, 
in particular 

^ {f ix)^{y)^{z)) - eSix - z)g{y, z) - e5{x - y)g{y, z). (20) 

We can write the corresponding equation for the vertex function, using (|17p : 
d f 

-^-^ J d*ud'^vig{y,u)r^{x;u,v)ig{v,z) = S{x ~ z)g{y - z) ~ S{x - y)g{y ~ z). (21) 

This form is easy to rewrite in the two-component formalism, taking into account that to satisfy the delta function 
the time arguments must be on the same contour. One finds in Fourier space 

fc^r^,,(fc;p,g) = [6abg^,\q) - 6acg^,\p)] i2n)^6ik+p-q). (22) 

In the Bloch-Nordsieck model, because of the special property of the vertex function expressed in p8)) . the vertex 
function is completely determined by the fermion propagator: 



^abc{k;p,q) = [(5afc^ft/(g) - SacgbJ'ip)] 



(23) 

p—q—k 
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III. SOLUTION OF THE DYSON-SCHWINGER EQUATIONS 



Since the vertex function in the Bloch-Nordsieck model can be expressed with the fermion propagator, the Dyson- 
Schwinger equations for the fermion propagator become closed. At zero temperature it can be shown |6[ that the 
solution of this equation yields the same result as the functional techniques, moreover, renormalization can be fully 
controlled here. In this Section we discuss the solution at finite temperature. 

We will use Feynman gauge, and denote the photon propagator as G^^, = —g^pG. Then this closed equation can 
be written as 



a',b' = l-^ ^ " 
.a — 1 ' 



Sa'b'iG ^)b'c{p) 



where U 



-ie U aa 



u^. In particular 
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La' = l 



'{k)Qaa'{p - k) - Sat 



'^^^ \G,a'{k)g,a'{p-k)- 



Sa'ciG ^)b'c{p-k) 

(27r)4 uk "-"-^ ' 



d^k Gn{k) 



(24) 



(27r)4 uk 

U^i:{Q-')a'2{p) J ^^G,Ak)g,.{p~k) 



(27r)4 uk 



(25) 



Instead of Gn and G22 it is more aesthetic to work with the retarded and advanced propagators (the relations 
are given in Since in the R/A formalism Gaa — 0, the retarded propagator satisfies a homogeneous self-energy 
relation 



Graip) = G(°)(p) + G(°)(p)I],,(p)G™(p), 

while the propagators in the 1, 2 components mix. From the definitions we easily find 

Sar — Sii + S12, Gil — G12 — Gar- 

Therefore we have, using ([25]) and ([6]) 

^arip)^Jip)G-a\p)~^M, 



(26) 
(27) 
(28) 



where 



J{P) 

AM : 



(27r)4 uk 
d^k Gii(fc) 



{G2l{k)Gra{p - fc) - Gra{k)Gl2{p-k)) , 



(2^)^ 



(29) 



It is easy to see that AM = 0. The 11 photon propagator Gn is even for k —k, which is true in general, but 
now we can prove by inspecting the free propagator which is exact in our case 



*Gii(fc) 



+ (n(fco) + e(-fco))27r sgn(fco) S{k^). 



(30) 



For the first term the k — >■ —k symmetry is evident, in the second we should use the identity n(feo) + n(— fep) + 1 = 0. 
Therefore with the change A; — >■ — fc of the integration variable, the Gn propagator remains the same while uk changes 
sign, so AM changes sign, too. As a consequence AM = 0. 

The Bloch-Nordsieck model, as all 4D interacting quantum field theories, contains divergences. To obtain finite 
result, we need wave function, mass and coupling constant renormalization. Since the above expressions have contained 
the original parameters of the Lagrangian, we should rewrite them in terms of the renormalized quantities. From now 
on the parameters m and e will denote the renormalized ones, while mg and cq are the bare quantities. Renormalization 
goes like in the zero temperature case @: assuming that the renormalized mass m = Zttiq where Z is the fermion 
wave function renormalization constant (this is ensured by the Ward identities) we can write Gra = Z{up — m) — Tiar, 
and from (l28t we find 



(31) 
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A. Calculation of J 

In the expression of J in eq. (|29)) there appears Qnik). As we discussed earher, for the sake of physical and 
mathematical consistency of the model, we must assume that the fermion describes a hard probe, itself is not a 
dynamical field, which means that we must set ^?i2(fc) = 0. Then from (1291) we can easily recover the zero temperature 
result 0. At finite temperature we have 

J[y) = -^e'U'J ^ ^G2i(fc)g™(p- fc). (32) 

Next we prove by recursion that the solution for Qra depends solely on w — up — m. It is true at tree level where 
Sra — up — m. So let us assume that Graip) — Qra{up — rn). Then 

J{p)^-ie'U^ j ^^^G2i{k)gra{np-m-uk), (33) 

implying J^{p) — J (up — m). Equation pip tells us that if depends only on up ~ to, then Qra also depends only 
on up — TO. With this statement the recursion is closed. 

Since in the BN model the free photon propagator is exact, we shall write it into eq. ([5^ . Using Q for the 
propagator, and applying the Landau prescription (w — > w + ie) we find 

Jiw) = e^U^J ^ (1 + n(fco))^((5(fco - k) - <5(fco + k)) Graiw - uk). (34) 

This result, as we shall show in Section UlI El is consistent with the results of 0, [I]- 

The k integration can be performed, apart from the single component q = ku. We find after a straightforward 
calculation: 

J{w) ^ I dq f{q,u)gra{w - q), (35) 

where a = e^/(47r) and 

uo{l+v) 

uo{l-v^) f ds q u^{l~v^) , eW("o(i~^)) _ ^ 

uo{l~v) 

where u = wo(l,v) and v = |v| (i.e. v is the velocity v = u/uq). 
At zero temperature /(<?) = 6 (<?)■ At z; = we find 

fiq,u = 0) = l + n{q). (37) 

B. Renormalization 

In (pSj) we find ultraviolet (UV) divergences. From the expression of f{q,u) (eq. ([Ml)) we see that for large 
momenta the thermal distribution functions always decrease exponentially, thus yielding UV finite result. So all the 
UV singularity is in the T — part, discussed already in 0- 

To apply the renormalized treatment at finite temperature, we recall some results from 6]. At T = psp can be 
written in spectral representation and with dimensional regularization as 



/9 / \ I 1 ^ / \ ^ I dw' f , 1 a f dw' _. 

Jo{w) = dq GraKw-q) = - — g{w ) dq — ■ — = - / — g(w ) 

TT J TT J ItT J q + W — W — te TT J ZTT 

— oo — oo 



oo oo _ 

„ , w' — w — ie 
V, - In 



(38) 



where 



2?e = ^ + ^ln(47r) + ^A/2 (39) 
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{Pi/2 = —1.96351 is the value of the polygamma function with 0, 1/2 arguments). 

As we discussed in ^] , the divergent term is necessary for the couphng constant and wave function renormahzation. 
We can write, assuming normahzability of g 



i + j{w) _ ~4 



—J- + X>e + J fin (w) 



(40) 



where J^finiw) is finite. We introduce 



47r2 4^2 

— +Ve = —., 

4 e2 



(41) 



where Zr and e now are finite (renormahzed) values. Using renormahzation group invariance we can write for the 
complete finite temperature contribution 



<(«■> - 5? 



dw' _ , A 
— g{w ) In — 

27J- 11^' 



w' — w — le 



- i dq{f{q,u) -<d{q))gra{w - q) 



(42) 



where e is a RG invariant coupling, A = /iexp(^^) is the momentum scale of the Landau-pole. The derivative of the 
first term reads 



w' — w — le 



2iT w — w' + ie 



-Graiw). 



(43) 



The imaginary part of I{w) term is zero for w < 0, moreover for w = it is negative (at least for large A), while for 
w — )■ — oo it is positive. So there exists a value w — ~M for which it is zero. Then we can write: 



I{w) 



dq Graiq)- 



(44) 



-M 



The scale M replaces the scale A. Assuming that M ^ T we can change the integration limits to —M M in the 
second part, too. Then we find 



CM 



47r2 



dqf{q,u)g{w - q), 



(45) 



where the integral symbol means J = J^j^,j- If it does not cause problem, we will send M — > oo. The zero temperature 
part is the same as in our earlier publication [6]. 
Summarizing, the renormahzed equation reads now: 



wg{w) 



dqf{q,u)g{w - q), 



(46) 



where f{q,u) is given by p6|) . Since this equation is linear, the same will be true for the spectral function (with 
different normalization conditions) 



wg{w) 



dqfiq,u) g{w - q), 



(47) 



C. Zero velocity case 

For V = 0, uq = 1 we find for (gT]) 



wg(w) ^ j dq{l + n{q))g(w - q). (48) 

TT 
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By sending the hmits of the integration to infinity, we reahze that the right hand side is a convohition. Therefore 
we change to Fourier space where it becomes a product, and the left hand side will be idtg{t). Using the Fourier 
transform of 1 + n{q) 



dw 
2^ 



-iT 



2 tanh(7rtT) 



we obtain the differential equation 



This has the following solution: 



idtg{t) 



iToi 



tanh(7rtT) 



git) = go (sinhTT^T)"^'^ 



(49) 



(50) 



(51) 



Before we proceed, we shall discuss this result. First we can easily recover the T = result, since for t ^ ^ the 
sinh function can be approximated linearly, and we get g(t) ~ t"^'^ . On the other hand this result is rather weird, it 
describes forever increasing correlation instead the physically sensible loss of correlation. Since this happens also at 
zero temperature, this is not an artifact of the finite temperature calculation. In accordance with Blaizot and lancu 
0, [1] , we should not consider this expression as the physical response function. Mathematically we can argue that we 
are not in the physically sensible analytic domain, the time dependent spectral function is not square-integrable for a 
real a value, as it should. We must therefore go over to the physical analytic domain, where the Fourier-transformation 
is well defined. 

For the analytic continuation we think equation ([5T|) valid as long as it yields sensible formulae, which is the case 
of imaginary a values. With this assumption the spectral function in the Fourier space will be an analytic function in 
a. For real a values the spectral function will be interpreted as an analytic continuation. We will see that it indeed 
provides sensible results. 

To perform the inverse Fourier transformation we apply Laplace transformation. With s± = ±iw we find 



J dte''"'g{t)= j dte-'-'g[t) + j dte''+'g{-t) = g+[s-) + {-^If/'' g+{s+), 



(52) 



where 



oo 

g+{s) = j dte-'' (sinhTTiT)"/" 



1 + ^ 

TT 



2tt 



a 
2^ 



f3s a 



Since the F-function satisfies r(l — z)T(z) — 7r/sin7rz, we can write with s = ±iw: 

\ 27r 27r 



r 1 



27r 



a 
2^ 



s—:Liw 



\ 27r 27r 



a Pw 

sin I — =F « 

' 2 ^ 2 



Then we get for the spectral function 

g{w) 



cosh(/3w) — cos a 



1 



r 1 



a Bw 

h I — 

27r 2tt 



where Na is a numerically determined normalization constant required by the sum rule 



(53) 



(54) 



(55) 



(56) 



to be satisfied. On Fig. [5] we can see the shape of the spectral function for different a values and for different 
temperatures. 

To discuss this result we make the following observations: 
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w/T w 
a.) b.) 

FIG. 2. a.) The exact, normalized spectral function at ?; = 0. Common features are the dominantly exponential decrease for 
w — >■ — oo, power-law decrease ~ w~^~°'^^ for lu — >■ oo and at the peak a finite curvature ~ a. b.) Temperature dependence of 
the spectral function at u = 0. In the limit T — >■ it is singular at ui = point. 



• 0{w) is a function of f^w only, which is understandable, since there is no other scale in the system which could 
form a dimensionless combination. 



• For a — )■ 0, we find 

2Trd{w), (57) 



cosh(/3u>) — cos a 

so we recover the free case. It is interesting, that this behavior periodically returns for a = 27rn. 

• For large values of w which is equivalent to the small temperature case we can use the asymptotic form of the 
F function for complex arguments with large absolute value: 

T{x) = e'^'x^ (x-i/^ ^ 0(^-3/2)^ (53) 

Then we find, up to normalization factors 

g{Pw » 1) — — Q{w)w-'-f. (59) 

cosh[pw) 

This is the well-known exact solution by Bloch and Nordsieck at zero temperature. Note that the 8 function 
came out correctly from the formula. At finite but small temperatures, for negative arguments we observe 
exponential decrease. 

This form also shows how at zero temperature we obtain zero wave function renormalization factor. The 
normalization factor (c.f. ([55])) is proportional to /3, while the asymptotic form is {I3w)^^^~ . Then approaching 
zero temperature we obtain T^w~'^~^ , which means a renormalization factor vanishing as ~ for T — >■ 0. 

• Now let us consider the w — > limit, i.e. the vicinity of the mass shell. We can expand g into power series 

, 4g(0)CT2 

Q{w) = -rj r , (60) 

{w - CTf + ( ^ - 1 ) C^T^ + 0{w^) 



where 



G 2 1 — cos a TT^ ZTT 



and ^(a) is the digamma function. The maximum of this function is at CT, the width is CT\/ AC^^ — 1. Since, 
however, the function is not symmetric, these parameters cannot be interpreted as a thermal mass and thermal 
width. For that we need to examine the real time dependence. 
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• For the real time dependence we use the fact that, according to ([55]) . q{p) = l3fo{l3{po — mj), which means that 
g{t) — e~™*fo{Tt). Omitting the oscihating phase (i.e. if we consider the envelope of g{t)), we recover the 
Fomier transform of /q. 

The real time dependence obtained from the inverse Fourier transformation of (j55p differs from (|5ip . This is 
because we performed an analytic continuation to the physically sensible analytic domain. The numerical inverse 
Fourier transform of the normalized spectral function (and, because iGra{t) — Q{t)g{t), for t > this is also 
the real time dependence of the retarded Green's function) can be seen on Fig. [31 




FIG. 3. Time dependence of the (envelope of the) retarded Green's function (or, equivalently, the spectral function) for v = at 
a = 0.5 on a logarithmic y-scale. For small times we find 1 — A(Tt)°'^^ , corresponding to the zero temperature time dependence. 
For large times it turns into an exponential exp(— aTt) damping. 



At small times we expect to recover the zero temperature result. Indeed, we observe g{t) = (1 — A(Ti)"/^)e^™* 
asymptotic form (for a — 0.5 this is valid up to Tt < 0.4), the power law time dependence is characteristic to the 
zero temperature result. At t = the value of the spectral function is 1, this is because of normalization. Note 
however, that naively at zero temperature we would obtain ~ faf-jr^-imt ^jjjjg dependence, describing growth of 
correlation and violating the normalization condition. Interpreting the zero temperature result as T — )■ limit, 
we could cure this apparent inconsistency of the model. At strictly T = we get back the physically sensible 
oscillating solution g{t) — e~™*. 

For large times (for tT > 1) the time dependence is ~ e'"'^*', which agrees with 0, [!]■ Comparing it to (|5ip we 
see that instead of an exponential rise we found an exponential decay, but with the same coefficient. This can 
be understood by noting that if we have a pole at it; = wq in the momentum space, meaning e~™°* exponential 
time dependence, this pole is present in the spectral function in position lOg, too. The physical retarded Green's 
function can have poles in the lower half plane, therefore we find in our case only the wo = —iaT pole, giving 
exponential damping. 

For the justification of the analytic continuation we also used a different method. We expanded the i-dependent 
result (|5ip into power series using 



(sinhrr)"/" = " X] ("\'''-^-''^ e'^'' + e(-x)(-l) ' 



k=0 

Now the inverse Fourier transformation acts on a pure exponential function. We use the formula 



(62) 



1 



iw 



(63) 



which is true, of course, if s > 0, but this is the formula for the analytic continuation, too. 
Fourier transformation, with appropriate normalization to ensure reality of g: 



Then the result of the 



e(H-E(-i)'(r. 



fc=0 



Sfc + iw Sk — iw 



(64) 
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where Sk = nT{2k — a/n). Using the 1)"/^'^ — cos §■ + i sin ^ definition we find after a simple calculation 



g{w) 



fc=0 



(65) 



The sum converges fast, and we can compare the result of tlie two calculations on Fig. 2] We can see that the two 



0.1 



0.01 



0.001 



0.0001 




-2 -1 1 2 3 4 5 
w/T 



FIG. 4. Comparison of the logarithm spectral function at a = 0.5 calculated from eq. (|55p and from eq. (|65l) . The two results 
agree well. 

methods of analytic continuation yield consistent result in the central peak regime. To understand the small deviations 
at the edges, we remark that if a is real then g{t = 0) = (c.f. (|51l) ). In Fourier space this means / dwg{w) — 0, 
therefore it can not be positive for all momenta. Using the second method, the position where the spectral function 
turns into negative values is, fortunately, at large Iw/rj values, therefore the peak is unaffected. The only precursor 
of the sign changing is the slight decrease at the edges of the plot. In our first method we started from imaginary a 
values, where lim g{t) ~ go =/= 0, then the normalization does not require negative values for g{w). 



D. Finite velocity case 

If t; 7^ we find for eq. (|T71) using the following formula 



uo{l+v) 



wg{w) 



a uo{l — V ) f ds 



IT 2v 



/ dq{l + n{-))g{w-q). 



(66) 



uo{l-v) 

The right hand side is again a convolution, and formally we can use the same method as in the v — case. We find 



"o(l+'j) 

„ auo{l-v'^) f ds Ts _ 
dtgit) = / -J — -TT-^TTTTt'W 



TT 2w J tanh(7rirs) 

«o(l-«) 



(67) 



which has the solution 



g{t) = g{0) exp 



"o(l+i') 
a uo(l — f ^) f ds 

Uoil-v) 



ln(sinh ntTs) 



(68) 



We cannot perform analytically neither the integral, nor its Fourier transform. But we can determine many features 
by investigating the i — >■ and i — > oo limits. 
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1. The limit t^O 



Since lim sinh TrtTs = ■KtTs: 

t->o 



"o(l+f) uo{l+v) 

lim "°^'-^'^ / ^ln(sinh7r^r.) = "°^^-^'^ / ^ InvrtT. = InTrTt + const., (69) 

uo(l— ") uo(l— ") 

where the constant comes from the integral of In s. Being a finite quantity, it goes into the normahzation. After 
exponentiation we find 

Q{t) ~ (T^)«/^ (70) 

which is the zero temperature result. So, as we expected the short time or large frequency regime reproduces the zero 
temperature case, and thus it is velocity-independent. 

2. The limit t — >■ oo 

Here the sinh can be approximated by the exponential, and so 

lim In sinh izTts = wTts - In 2 (71) 

The In 2 yields a constant factor which goes into the normalization. The rest gives, including the prefactors 

uo(l+v) 

^ Itnfl — W^) f dS , sr^ 

aTt \^ ' J -=a,jf{u)Tt, (72) 

uo(l— v) 

where 

aeffiu) = a In J . (73) 

From this form we obtain for the spectral function in the asymptotic limit: 

^(t) = Ce«=//(")^*. (74) 

We can easily check that lim ae//(w) = a. Therefore the u — >■ limit is analytic. 

Since in the asymptotic time regime we simply get the substitution rule a — > aeff{u) as compared to the v = 
case, the analysis of the vicinity of the peak of the spectral function, and the large time dependence will remain valid 
in the finite velocity case, too. with a modified value of the coupling. In particular, since «£//(«) < a, we obtain a 
smaller damping, larger lifetime for v > cases. Physically this property is the consequence of the decreasing cross 
section at larger energies. 

3. Solution for t € (0, oo) 

For intermediate times we could not work out analytically the integral. Nevertheless, we have a well-controlled 
numerical method to find the spectral function, once the analytic behavior for large t is identified. We express the 
wanted Quit) as a product of the known Qu=o{t; o^ef f) and a correction factor 

Z{t)e^=o{t;aeff), (75) 

where 



exp 

Z{t) ^ 



uo{l+v) 

auo{l — v"^) f ds 



2lTV J 

uq(1-v) 



ln(sinh7rtTs) 



(sinhTrT^)^^ 



(76) 
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After a short algebra we find 

2^ 7 (sinh nTt)' [ ' ^ ^ 

uo{l-v) J 

o-i»eff (") 

The so-defined ratio is symmetric Z{t) = Z{—t). For small t arguments it behaves as Z{t) ^ (Tt) s and, since 
a ^ aeff{u), we also know Z{t = 0) = 0. At large t we find lim Z{t) = 1. We can determine it numerically, for a 

specific V it can be seen on Fig. [5] 



Z{t) = exp ■ 



1 



10 



100 



0.01 



0.1 



FIG. 5. The Z{t) function on a logarithmic plot. For small times it is a power, for larger times it flattens out. 



We can numerically Fourier transform Z(t), and perform a convolution in the Fourier space with the gy=o{w) 
function (jSTj) . This ensures that we use the same analytic continuation for the different velocity cases. As a result we 
obtain Fig. |6l We can observe that the peak becomes narrower for larger velocities, corresponding to the decreasing 
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0.01 



0.001 




r|=0.0 
Ti=0.96 
r|=1.31 
r|=2.06 
ri=2.77 



a.) 



b.) 



FIG. 6. Velocity dependence of the spectral function for a = 0.5. The 77 values are rapidities, v = tanhij. a.) is a linear-linear 
plot to show that the peak region becomes more and more peaked with increasing 77 (u). The log- log plot demonstrates that 
the asymptotics remain the same. 



aeff value. At large momentum the asymptotics is the same for all velocities (for a given a), because the zero 
temperature result is insensitive to the value of v. 

Here again we can work out the real time dependence. We now write Quip) = uo/3/u(/3(poWo — pu — to)) and find 



fT 

Mo 



(78) 



The result of the numerical inverse Fourier transform can be seen on Fig. [T] At small times the spectral function 
(retarded Green's function) is velocity independent, this is the zero temperature asymptotics. For large times (for 
Tt > 1) the time dependence turns into ^ Q-o'„ff{u)Tt_ 
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FIG. 7. Comparison of the real time dependence of the retarded Green's function (or, equivalently, the spectral function) 
for zero and finite velocity at a — 0.5 on logarithmic y scale. The small time behavior does not change, at large times the 
exponential damping turns to exp{—aeff{u)Tt). 



E. Discussion of earlier results 



We can compare our results to the earlier results in the literature. The long time asymptotics of the finite tem- 
perature solution of the Bloch-Nordsieck model was already discussed in 0, They followed a different, functional 
approach. Still, the two methods lead to the same intermediate result. Neglecting renormalization effects (which is 
treated later in Ref. 0, @|), our eq. ([31]) together with ([M)) yields, using the notation w — up — m 



wQra{w) = e^U^ 



d'^k 1 



(27r)4 uk 



(1 + n{ko))g{k) Graiw - uk), 



(79) 



where g^k) — ff (^(fco ~ k) — S{ko + k)) is the photon spectral function. After Fourier-transformation we find expo- 
nentiation of the real time contributions 



Grait) ^ Grait = 0)e^'^'\ 



(80) 



where 



AJL(^l + n{ko))g{k)- 



-iukt 



[uky 



-it$o + + lnA(i), 



where $0i ^{t) and InA(t) are real quantities, corresponding to the notation of [8]. Using = \ this means 



m 

A 



d'^k g{k) 



ikt 



ikt 



d-^k g{k) 



(27r)4 uk 



exp 



-^^4 n{ka)g{k) — 



- cos uki 



(82) 



These expressions agree with the equations (2.24) and (2.25) of (the constant phase $o has no physical meaning). 

The analysis of this formula, however, differs in our case and in vl, Isj- We strictly restrict ourselves to the original 
Bloch-Nordsieck model, and used the free photon spectral function. In Ref. 0, [1] the authors used HTL-improved 
photon spectral function (cf. their eq. (3.1) and (3.2)). As it turns out, the most important contribution comes from 
the small frequency limit of the continuum (Landau damping) part. This explains why the asymptotic time behavior 
differs in our case and in the case of Ref. 0[1] ( exp(— ae//(w)rt) vs. exp(— Cilogt) ). 

In Ref. [lo| Fried et. al. use again a different formalism. Since they examine a different physical situation, the 
comparison is much more difficult. What is clear, however, that they also use the original version of the model, and 
also find exponentially dampin g so lution. 

It is very interesting that in IJj the authors found the same exp(— aTi logi) like solution as was the case in Ref. 
0, Q, although with a different line of thought. They use the dynamical renormalization group idea fis'l, where the 
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secular terms are melted into finite time dependence of the renormalized parameters. Clearly they cannot consider 
all photonic diagrams, just those which contribute to the renormalization group (RG) equations. The logarithmic 
enhancement of the damping there can be interpreted physically as an eternally growing cross section of the incoming 
hard particle which collects more and more soft photons around itself. The analysis of the pure Bloch-Nordsieck 
model results in a finite damping, which means that in this model the initial growth of the cross section eventually 
stops, the soft photon cloud saturates. The physical interpretation of the saturation probably is that the multi-photon 
contributions arriving from different spacetime points become incoherent. 

The two scenarios, one with ever growing photon cloud, the other with saturation, are both approximations of the 
real QED (RG, HTL and free photon approximations, respectively). The question that which one is finally manifested 
in QED, can be answered only after a full analysis of the complete QED where all these effects are present. 



IV. CONCLUSIONS 

In this paper we studied the Bloch-Nordsieck model at finite temperature, in particular we studied the fermionic 
spectral function. We used the strategy introduced in 6] which is based on the Dyson-Schwinger equations, where 
the infinite hierarchy is closed by using the Ward identities for the vertex function. We worked out the corresponding 
equations at finite temperature in the real time formalism, and solved them. This procedure is exact in the Bloch- 
Nordsieck model. 

At zero velocity we were able to obtain fully analytic results for the spectral function. For large momenta and/or 
zero temperature this formula agrees with the zero temperature result. At finite temperature there appears an 
asymmetric peak which decreases exponentially below the mass shell {u^p^ < m) and as a power law above the mass 
shell (u^p^ > m). 

We also worked out the real time dependence which has two characteristic regime. For small times, starting 
from its initial value, it behaves as a power law Q{t) = (1 - A(Tt)"/'')e-*"*, where A depends on a. The naive zero 
temperature calculation yields ~ fo/Tr ^-imt ^jj^^g dependence which is not normalizable and corresponds to a physically 
hardly interpretable forever growing retarded response function. With the finite temperature as a regulator, we could 
interpret the zero temperature result, and we got a physically sensible purely oscillating response function. 

For large times we find exponential damping. The damping rate is aeff{u)T, where the effective coupling aeff{u) 
is given in (I73p . The damping is smaller, the lifetime is longer for larger velocities, which physically can be interpreted 
as the consequence of decreasing cross sections. We remark that the damping in the pure Bloch-Nordsieck model 
differs from the one with HTL-improved photon propagator; in this latter case one finds a faster-than-exponential 
damping with an exponent ^ — tlni. 

We expect that the method we worked out for the Bloch-Nordsieck model can be applied, as an approximation 
scheme, also for the full QED. Hopefully the renormalizability of the resummation will remain true, too. 
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